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Abstract — Convolutional network-error correcting codes 
(CNECCs) are known to provide error correcting capability 
in acyclic instantaneous networks within the network coding 
paradigm under small field size conditions. In this work, we 
investigate the performance of CNECCs under the error model 
of the network where the edges are assumed to be statistically 
independent binary symmetric channels, each with the same 
probability of error p e (0 < p e < 0.5). We obtain bounds on the 
performance of such CNECCs based on a modified generating 
function (the transfer function) of the CNECCs. For a given 
network, we derive a mathematical condition on how small p e 
should be so that only single edge network-errors need to be 
accounted for, thus reducing the complexity of evaluating the 
probability of error of any CNECC. Simulations indicate that 
convolutional codes are required to possess different properties 
to achieve good performance in low p e and high p e regimes. 
For the low p e regime, convolutional codes with good distance 
properties show good performance. For the high p e regime, 
convolutional codes that have a good slope (the minimum 
normalized cycle weight) are seen to be good. We derive a lower 
bound on the slope of any rate b/c convolutional code with a 
certain degree. 



I. Introduction 

Network coding as a means of increasing throughput in 
networks has been extensively studied in IU-0. Block 
network-error correction for coherent network codes has been 
studied in H-16). In all of these, the sufficient field size 
requirement for designing good block network-error correcting 
codes (BNECCs) is quite high. To be precise, the sufficient 
field size requirement for constructing a BNECC along with a 
network code which corrects network-errors due to any t edges 
of the network being in error once in every J network uses is 

such that q > |T| J '" 
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where T is the set of sinks. This 



requires every network-coding node of the network to perform 
multiplications of large degree polynomials over the base field 
each time it has to transmit, and therefore is computationally 
demanding. Moreover, the bound increases with the size of 
the network. It is therefore necessary to study network-error 
correcting codes which work under small field size conditions. 

Convolutional network-error correcting codes (CNECCs) 
were introduced in [7] in the context of coherent network 
coding for acyclic instantaneous networks. The field size 
requirement for the CNECCs of [7| is independent of the 
number of edges in the network and in general much smaller 
than what is demanded by BNECCs. Although the error 
correcting capability might not be comparable to that offered 
by BNECCs, the reduction in field size is a considerable 
advantage in terms of the computation to be performed at each 



coding node of the network. Also, the use of convolutional 
codes permits decoding using the Viterbi decoder, which 
is readily available. CNECCs with similar advantages for 
memory-free unit-delay acyclic networks were discussed in [8] 
and the benefit obtained in the performance of such CNECCs 
by using memory at the nodes of unit-delay networks was 
discussed in |9j . 

The CNECCs of [7] were designed to correct network-errors 
which correspond to a set <f> of error patterns (subsets of the 
edge set) once in a certain number of network uses (a network 
use being the use of the edges of the network to transmit a 
number of symbols equal to the network code dimension). 
A similar error model (with $ being all subsets of the edge 
set with t edges) was considered in |4|-|6|. While this error 
model allows code construction, it is less realistic because the 
errors corresponding to any error pattern in $ are assumed to 
occur with equal probabilities. 

A more realistic error model would be to assume every 
edge e in the network as a BSC with a certain cross-over 
probability (p e ) and with errors across different edges to be 
i.i.d. In this paper, we assume such an error model (with 
p e being the same for all edges) and analyze CNECCs over 
the binary field. Binary network codes together with this 
error model were studied in 1 101. The decoding of BNECCs 
under a similar probabilistic setting was discussed in IfTTll . 
However, practical analysis and simulations of BNECCs under 
a probabilistic error setting is difficult because of the large field 
size demanded. On the other hand, the CNECCs developed in 
require small field sizes and thus facilitate analysis. The 
contributions and organization of this paper are as follows. 




Fig. 1. Butterfly network 

• After briefly discussing CNECCs for the network coding 
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setup (Section |nli, we present the error model for the net- 
work. If the edge cross-over probability p e << 0.5, then 
it is sufficient to compute only single edge network-error 
probabilities in the network thereby reducing the com- 
putations required to study the performance of CNECCs. 
For any network with a given number of edges, we derive 
a bound on how small this p e should be so that this 
assumption of ignoring multiple edge network-errors can 
be made safely. (Section [TTTb 

Expressions for the upper bound on the bit error probabil- 
ity of CNECCs are obtained based on a modified version 
of the augmented path generating function (T(D,I)) of 
the CNECC being used. (Section HVb 
We analyze the performance of CNECCs on networks 
with a probabilistic error model using simulations with 
the butterfly network (Fig. Q]) as an example. Simulations 
on the butterfly network indicate that different criteria 
apply for CNECCs to be good under low and high 
p e conditions. We therefore suggest different types of 
CNECCs under these two conditions. (Section [Vli 
For high p e conditions, it is seen that those codes perform 
better which have a high value of slope, which is defined 
as follows. 

Definition 1 ( 4721/ ): Given a minimal encoder of a rate 
R = b/c convolutional code C, the minimum normalized 
cycle weight 



w„(o) 
n := mm < H \ 
oEO\ 0l { l(o) 



(1) 



among all cycles o G C(the set of all cycles) in the state 
transition diagram of the encoder, except the zero cycle o\ 
in the zero state, is called the slope a of the convolutional 
code C. Here w H (o) indicates the Hamming weight 
accumulated by the output sequence while traversing the 
cycle o, and 1(6) is the length of the cycle in c-tuples. 
• We derive a lower bound on the slope of any rate b/c 
convolutional code over any finite field (Section [VB , and 
conclude with a short discussion of the paper and several 
directions for future research (Section [Vill i. 

While CNECCs only over F2 are considered for the analyses 
and simulations of this paper, CNECCs over any field size can 
be studied using similar methods. 

II. Convolutional codes for network-error 

CORRECTION 
A. Network model and network code 

An acyclic network can be represented as an acyclic directed 
multi-graph (a graph that can have parallel edges between 
nodes) Q = (V, £) where V is the set of all vertices and £ is 
the set of all edges in the network. Every edge in the directed 
multi-graph representing the network has unit capacity (can 
carry utmost one symbol from F2). 

Let n be the mincut between the source s and the set 
of sinks T and the dimension of the network code. An 
n-dimensional binary network code can be described by 
three matrices A (of size n x \£\), F (of size \£\ x \£ |),and 



B T (of size \£\ x n for every sink T G T), each having el- 
ements from F 2 . Further details on the structure of these 
matrices can be found in 0. The network transfer matrix 
corresponding to a sink T is an n x n binary matrix Mt 
such that for any input x G F 2 , the output at sink T G T is 
xM T = xAFB T . 

B. CNECCs 

For a given set of error patterns $ and for some k < n, 
a method of constructing rate k/n convolutional codes was 
given in Q such that these CNECCs will correct network- 
errors which correspond to the patterns in $. For a given 
network with a network code, the definitions for the input 
and output convolutional code are as follows. 

Definition 2: An input convolutional code, C s , correspond- 
ing to an acyclic network is a convolutional code of rate 
k/n(k < n) with a input generator matrix Gj(z) implemented 
at the source of the network. 

Definition 3: The output convolutional code Ct corre- 
sponding to a sink node T in the acyclic network is the k/n 
convolutional code generated by the output generator matrix 
Go,t{z) which is given by Go,t{z) = Gj(z)Mt, with Mt 
being the full rank network transfer matrix corresponding to 
an n-dimensional network code. 

It was shown in Q that errors corresponding to $ can be 
corrected at all sinks as long as they are separated by a certain 
number of network uses. Moreover, a sink can achieve this 
error correcting capability by choosing to decode on either 
the input or the output convolutional codes depending upon 
their distance properties. 

Example 1: Table U shows the network transfer matrices of 
the butterfly network of Fig. Q] and an example of a CNECC 
along with the output convolutional codes at the two sinks. 

TABLE I 

Butterfly network of Fig.Q]with the input convolutional 

CODE G/(z) = [1 + z + z 2 1 + z 2 ]. 



Sink 


Network transfer 
matrix 


Output convolutional code 


Ti 


M Tl = 


[o l) 




Go,T 1 (z) = [l + z + z 2 z] 


T 2 


M T2 = 


(! 




Go,T 2 (z) = [z 1 + z 2 } 



III. Network-errors in the BSC edge error model 

Any edge e G £ in the network is assumed to be a binary 
symmetric channel with probability of error being p e and 

errors on different edges are assumed to be i.i.d. A network- 

\£\ 

error is a vector w G F 2 with Is at those positions where the 
corresponding edge is in error. The probability of a network- 
error KJgFf is ±en P e H{w \l -p e )\ £ \- w H( w ). 

Let e T denote the random error vector at sink T. The 
probability that e T = y G F 2 l is as follows. 

o H {w) 



(y) 



E 



Pe 



'(l-p^l^l-^M (2) 



\£\ 



:wFB T -- 



yPl(l-p e ) 



£\- 
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where a; „ indicates the number of network-error vectors from 

lei 

¥ 2 with weight i, such that they result in the error vector y 
at sink T. 

For any given network, it is essential to calculate the error 
probability of e T being any y £ for each sink T € T 
in order to analyze the performance of any CNECC over 
the network. Equation (f2]) indicates that this involves a large 
number of computations even if the given network is small. 
However, if p e << 0.5, then it is sufficient to compute 
only single edge network-error probabilities for any particular 
error vector at any sink, thereby reducing the number of 
computations. In particular, suppose 



a 1>y p e (l -p e ) |£:| 1 > A VV y p*(l -p e ) 



\£\ 



(3) 



for any error y at any sink T with a\ y ^ 0, for some A > 0. 
We then have the following upper bound. 

(4) 

The probability of the error vector being £ F£ is upper 
bounded independent of A as follows. 



Pe (0) < 1 



E 

»eB?\{0}:oi, w5 40 



ai, y Pe(l - Pe) |£|_1 (5) 



If p e is small enough so that (01 holds for some large A, then 
the upper bounds of (HJ and (0 become tight, and hence single 
edge network-errors alone can be considered in the network 
without any significant loss of generality. 



A. An upper bound on p e 

In this subsection, we obtain a sufficient upper bound on 
p e for a given network for (0 to hold so that only single 
edge network-error probabilities need to be calculated. This 
bound obtained holds for any network with a given number 
of edges and is independent of the network code chosen. It 
is seen that this bound on p e is inversely proportional to the 
number of edges in the network. This is a reasonable result 
because among the network-errors which result in some error 
vector at a sink, the difference between the number of multiple 
edge network-errors and the number of single edge network- 
errors would in general increase with the increase in network 
size, thus lowering the value of p e upto which 01 would hold. 
Towards calculating this bound, we first prove the following 
lemma. 

Lemma 1: For any integer m > 1 and V0<p< 1, 

(1 - p) m >l-mp. 

Proof: For any p > ~, 1 — rap < 0, and the proof is 
obvious. Therefore we prove the lemma only for p < ^. 
We have 



(1 - p) m -i-mp = Y j 



(-l)V 



(6) 



Let S be defined as 

E 



Therefore the R.H.S of © becomes 



E 



(-i)V 



S if m is odd 
S + p m if m is even 



If S > 0, the lemma is proved. Now every element inside the 
summation of S is of the form 



i+i 



(7) 



Since p < —.we have 



PM1-?E?P >o, 



in — z \ / m — i 1 , 
l-—^p) > [1-—^.- | >0. 



i + 1 



i + 1 m 



This means that every element in the summation of S is non- 
negative, which means that S > 0, hence proving the lemma. 

■ 

We now state and prove Proposition Q] which gives the upper 
bound on p e for (f3]l to hold. 

Proposition 1: For any error y at any sink T with ai, y ^ 0, 
the following holds 



\£\ 



ai, y pe(l - Pe) 181 1 > A ^ai )S p*(l -p e ) 



S\-i 



i=2 



if 



Pe < 



(|£|-1)(A|£|-A + 1)' 



Proof: Since a,i, y < 



\£\ 



\/i and a\, y > 1 corre- 



sponding to such an error y as considered in the proposition, 
it is sufficient to consider the following case 

Pe (i- Pe f^>\(f:( f )pi(i-p e y £ ^ 

to get the bound on p e . Hence, we have 

Pe (l -Pef^ 1 > A (l " \£\Pe{l -Pe)' £hl - (1 - Pef 1 ) 

(8) 

^(l-p e )l £ l- 1 ((A|f| + l)p e + A(l- Pe )) > A (9) 
By Lemma Q] the inequality of (0 holds if the following holds 

(1 - (\£\ - l) Pe ) {{\\E\ + l) Pe + A(l - Pe )) > A 

(1 - (\£\ - 1) p e ) (X\£\p e - (A - l) Pe + A) > A (10) 
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Probability of error vector [1 0] at Sink T 
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1 0*(Sum of probability of multiple edge 
network-errors that result in vector [1 0] at Sink T ) 

Probability of single edge network-errors 
" that result in vector [1 0] at Sink T 
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Computations involving 
multiple edge network-errors 
can be ignored if p g < 0.0135 




0.025 



0.015 



Probability of an edge being in error (p ) 



0.005 



Fig. 2. Threshold p e for the butterfly network at Sink Tj 



Simplifying dTOb . we get 

Pe < 



1 



(\£\-l)(\\£\-X + l) 

■ 

The bound of Proposition Q] holds for any network with 
\£\ edges for a chosen A and in general is loose as indicated 
by Fig. |2] Having chosen A = 10, Fig. [2] shows the single 
edge network-error probabilities and 10 times the multiple 
edge network-error probabilities obtained using simulations 
with respect to varying p e , corresponding to the error vector 
[1 0] at Sink T\ of the butterfly network. The threshold p e is 
approximately 0.0135, which is the lowest computed for any 
error vector at Sink T\. A similar value can be computed for 
Sink T2 . This is approximately an order of magnitude greater 
than what the bound of Proposition Q] indicates (p e < 0.00154 
for the butterfly network which has 9 edges). 

IV. Bound on the bit error probability of a 
CNECC 

We can bound the bit error probability of a CNECC follow- 
ing lTT3l upon a slight modification of its augmented generating 
function T(D,I), which is a polynomial in D and / where 
any element of T(D,I), say bD d P, indicates b number of 
paths which are unmerged with the all-zero codeword with 
a Hamming distance of d and i number of input Is being 
encoded into the unmerged codeword segment. We compare 
the bound thus obtained with simulations on the butterfly 
network in Subsection IV-BI 

However, because the network coding channel has F™ inputs 
and F™ outputs, the generating function of the convolutional 



code needs to be modified to capture every n bits transmitted 
at once. 

Therefore, we use the place-holders D v for the branches 
of the state transition diagram with the output vector being 
v G F™\{0}. The modified augmented generating function, 
T(Dqoo..qi, •••) -Din. .11, 1), is thus the transfer function of the 
convolutional encoder with the state transition diagram with 
the branches weighted with appropriate D V P. 

The bit error probability for a given rate k/n CNECC for 
a sink T is then bounded as 

1 dT(Dooo..oi, -Dili.. iij I) 
k dl 



< 



(ID 



1=1, D„ 



where 



Zv,T 



E 



Pe^ [y)p err (y + v) 



is the Bhattacharyya bound on the pairwise error probability 
between and v, with p &T (y) being the probability that the 
error vector obtained at sink T after applying the inverse of the 
network transfer matrix (My) is y. The partial derivative of 
( fTTT l can be upper bounded according to the numerical upper 
bound ([T2l i shown at the top of the next page. 

Example 2: Fig. |4] shows the state transition diagram cor- 
responding to a minimal encoder (controller canonical form) 
of the convolutional code generated by the matrix [1 + z + 
z 2 1 + z 2 ) . The modified augmented generating function can 



be obtained as 



T(D 01 ,D 10 ,D U ,I) 



ID 2 n D 10 



P (D 2 



in 



n 2 ) 
^01) 



21D 



10 



(13) 



5 



BER at SinkT 2 

0.5 




Probability of an edge being in error (p ) 



Fig. 3. BER at Sink T 2 for two CNECCs (Table Hp 



dT{D Vl ,...,D V2n _ t ,I) 



01 



I=l,D„.=Z v 



T(Z Vl , Z v , 1 + e) - T(Z Vl , Z v 1) 

< where e << 1. (12) 

e 




Fig. 4. State diagram of the code generated by [1 + z + z 2 1 + z 2 ]. 

It can be noted that with D v = D w h( v ) in (Qj]), the usual 
augmented generating function T(D,I) of the code can be 
obtained. 

V. Inference via simulation results 

A. Decoding of CNECCs 

Given a p e value at which the network operates, any sink 
can choose to decode a CNECC either on the trellis of the 
input convolutional code or that of its output convolutional 



code, depending on their performance at the given p e value. 
Decoding on the output convolutional code is advantageous 
to any sink because it does not have to perform the network 
transfer matrix inversion before having to decode every time 
it receives the incoming symbols. 

TABLE II 

2 CNECCS FOR THE BUTTERFLY NETWORK (FlG.Q} WITH THE OUTPUT 
CONV. CODES AT THE SINK T 2 . 



Input convolutional code 
generator matrix 


Output convolutional code 
generator matrix at Sink T2 


[1 + z 1] (C) 


[2 1] (Co,T 2 ) 


[i A (C) 


[1 + z z\ (C' O T J 



Example 3: Fig. [3] shows the performance of two CNECCs 
and their respective output convolutional codes (shown in 
Table Hill at sink T2 of the butterfly network. It can be noted 
that for all p e values shown, code C To performs better than 
code C . Thus if the code C is used, sink T2 can always decode 
on the trellis of C' T ^ . The opposite situation is observed for 
the pair C and Co.t 2 - It is therefore more beneficial for sink 
T2 to decode on the trellis of C (after matrix inversion) for 
any p e < 0.25. For p e > 0.25, sink T2 can decode on the 
trellis of Co,t 2 > as the performance improvement obtained by 
decoding on C is negligible. 
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cc 
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0.15 
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0.05 



Code C 1 (Free distance = 3, Slope = 1) 
Code C 2 (Free distance = 5, Slope = 1/2) 




P e = 0.025 (where the performance of the codes gets interchanged) 



5.3 



0.25 0.2 0.15 0.1 0.05 

Probability of error of an edge being in error (p ) 



Fig. 5. BER and bounds on BER at Sink T\ for 2 codes 



B. Coding for different values of p e 

Fig. [5] shows the performance of two different CNECCs 
(shown along with their properties in Table [TlTb at Sink T\ 
of the butterfly network. Similar performances are seen at 
Sink T2- The decoding for all these CNECCs are done on the 
corresponding input convolutional code. It is seen that there 

TABLE III 
CNECCS FOR THE BUTTERFLY NETWORK 



CNECC generator matrix 


Free distance 


Slope 


Gi{z) = [l + z 1J (Ci) 


3 


1 


Gi(z) = [l + z + z z l + z J ] (C 2 ) 


5 


1/2 



are two regimes of operation (for each pair of convolutional 
codes) where the performance of the codes get interchanged. 
This was already noticed in lfl2ll in the context of AWGN 
channels. The value of p e for which these regimes becomes 
separated is not only dependent on the CNECC-pair chosen, 
but also on the network and the network code, and would 
probably decrease with the increase in the size of the network. 

1) Coding for the low p e regime: Fig. [6] shows the per- 
formance of convolutional codes with different free distances 
on the butterfly network for low values of p e , along with 
the bounds on the bit-error probability evaluated according to 
Section [IV] Codes with better distance spectra are good in the 
low p e regime. According to Fig. [5] this behavior is seen upto 
p e = 0.025, however the bounds on the bit-error probability 
states become very loose beyond p e — 0.005 which is why 
the p e has been restricted to that value in Fig. [6] 

Maximum Distance Separable (MDS) convolutional codes 
thus seem to be a good choice. The design of such con- 



volutional codes along with the bounds on the field size 
requirement was discussed in Q for a fixed set of error 
patterns. If the value of p e is low enough, one might follow 
the design given in Q assuming the set of errors to be all 
possible single or double edge network-errors alone. 

2) Coding for the high p e regime: From Fig. [5] it is 
seen that codes with higher slopes are good for the high 
p e regime. The definition of the slope a of a convolutional 
code C is as in (HJ. For a given memory m and free distance 
dfree, a convolutional code is said to be a maximum slope 
convolutional code 1121 if there exists no other code with a 
higher slope for the same memory and same free distance. 
Families of convolutional maximum slope convolutional codes 
were reported in fl2l . discovered using computer search. 

VI. A LOWER BOUND ON THE SLOPE OF RATE b/c 
CONVOLUTIONAL CODES 

As seen in Subsection IV-BI codes with good slopes perform 
well in high p e conditions. It is therefore important to inves- 
tigate the properties of the slope parameter and to come up 
with constructions which yield codes with good slopes. Upper 
bounds on the slope of convolutional codes were given in IPT21 . 
lTT5l . A lower bound on the slope of any rate 1/c convolutional 
code was given in 11151 . In this section, we derive a lower 
bound on the slope of any rate b/c convolutional code over 
any finite field. 

A primer on the basics of convolutional codes can be found 
in Appendix [A] Towards obtaining a bound on the slope a, 
we first give the following lemma. The proof of the following 
lemma is on the lines of Lemma 1 in Q- 
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Fig. 6. BER and bounds on BER at Sink T\ in the low p e regime 



Lemma 2: Let C be a rate b/c convolutional code with 
degree 8. For some i > 0, if there exists a 6+1 length partial 
codeword sequence 



V[i,i+6] ■= [Vi,Vi+l, ...,V i+ s\ 



where Vj = £ for j = i, i + 1 , 



(5, then t> [lji+(5 ] has 



at least one cycle around the zero state of the corresponding 
minimal encoder of C. 

Proof: Let G m b(z) be a minimal basic generator matrix of 
C) . Let the ordered Forney indices (row degrees of G m t,(z)) 
be vi,V2, ■■■,i / b — v max , and therefore 8 being the sum of 
these indices. Then a systematic generator matrix(G SJ , s (z)) 
for C that is equivalent to G m b{z) is of the form 

G sys {z) = T~ 1 (z)G mb (z) 

where T(z) is a full rank 6x6 submatrix of G m b{z) with a 
delay-free determinant. We have the following observation. 

Observation 1: The degree of det (T(z)) is utmost 6. Also, 
we have the (i,j) th element U t j(z) of T^ 1 (z) as 



tij(z) 



Co factor (T(z) 



3^1 



det (T(z)) 



where Co f actor (T(z)j.i) 6 ^q[z] is the cofactor of the 
(j,i) th element of T(z). The degree of Co f actor (T(z) jti ) 
is utmost 8 — Vj < 6 — v\ . 

Let aij(z) £ ^q(z) represent the (i,j) th element of 



G sys (z), where 



ai,j(z) 



b 

E 

fc=i 



h,k{z)gk,j{z) 



_ J2k=i Cofactor(T(z) k ,i)gk,j(z) 
det(T(z)) 

9k,j(z) being (k,j) th element of G m b{z). Therefore, the 
element aij(z) can be expressed as 

Pi,j( z ) 



ai,j(z) 



det (T{z)) 



where the degree of Pij(z) e ¥ q [z] is utmost 5 + v rnax — v\. 
Now if we divide Pij(z) by det (T(z)), we have 



a t ,j{z) = q itj (z) + 



det (T(z)) 



where the degree of qij(z) £ F„[z] is utmost 



(14) 

, — v\, and 

the degree of r^iz) is utmost 8 — 1. Because every element 
of G sys (z) can be reduced to the form in ( fT4b . we can have 
a realization of G sys (z) with utmost 8 memory elements for 
each of the 6 inputs. Let this encoder realization be known as 

K 

^ sys • 

Now we shall prove the lemma by contradiction. Let v(z) 
be a codeword which contains the partial codeword sequence 

v [i,i+8] as follows: 

v(z) = [v ,v 1 , ...,Vi = 0,0, ...,0,v i+s =0,v i+ s+i,...] 

Let u s (z) be the information sequence which when encoded 
into v(z) by the systematic encoder E sys . Because of the 



x 




If the Hamming weight of the path 
in all the consecutive 5+1 segments 
is zero, then a = a =0. 

1+5 1+5+1 



Fig. 7. The trellis corresponding to a systematic encoder of C 



systematic property of E sys , we must have that 

U s ,i = w s ,i+i = •■• = u s ,i+s = £ Fj- 

By Observation Q] E sys is an encoder which has utmost 5 
memory elements (for each input), and hence the state vector 
(T i+ s e at time instant i + S becomes zero as a result of 5 
zero input vectors. Fig. [7] shows the scenario we consider. 

With another zero at time instant i + S, there is a zero cycle. 
But we need to prove it for a minimal encoder, not a systematic 
one. So, we consider the codeword v (z), which can now be 
written as a unique sum of two code words v(z) = v'(z) + 
v"(z), where 

i+S 

v'{z) = ^2 v k zk = bo, Vi = 0, v i+s = 0, 0, ...] 

fc=0 

and 

v"{z) = Vkzk = [0,0, ...,0,0, »(+*+!,...] 

k=i+8+l 

where OgFJ and the uniqueness of the decomposition holds 
with respect to the positions of the zeros indicated in the two 
code words v'(z) and v"{z). 

Let u m b(z) be the information sequence which is encoded 
into v(z) by a minimal realization E mb of a minimal basic 
generator matrix G m b{z) (a minimal encoder). Then we have 

u mb (z) = u' mb (z) + u'^ b {z) 

where u' mb (z) and u" nb (z) are encoded by E mb into v'(z) 
and v"(z) respectively. 

By the predictable degree property (PDP) [14] of minimal 
basic generator matrices, we have that for any polynomial code 
sequence v(z), 

deg (v(z)) = max {deg (u mb j(z)) + v{\ . 



where u mb j(z) S ^q[z] represents the information sequence 
corresponding to the I th input, and deg indicates the degree 
of the polynomial. Therefore, by the PDP property, we have 
that deg (u' mb (z)) < i, since deg (v'(z)) < i. 

Also, it is known that in the trellis of corresponding to a 
minimal realization of a minimal-basic generator matrix, there 
exists no non-trivial transition from the all-zero state to a non- 
zero state that produces a zero output. Therefore we have 
deg (tij^jz)) > i + 6 + 1, with equality being satisfied if 
v i+S+i 0- Therefore, u mb (z) is of the form 

u mb (z) = u' mb (z) + u" nb (z) 

2 — 1 oo 

= u 'mb,k zk + u mb,k zk 
k=l k=i+S+l 

u mb (z) = [u' mbt0 , .., u' mb ^j_x, 0, 0, ..] 

+ [0, .., 0, u mb i+s+1 ,u mb i+s+2 , ..] 

i.e, if 

U mb (z) — [u mb fi, U m b,l, —, U mb ,i, ...,U mb ,i+8,Umb,i+S+l, ••] 

then it mM = u mbii+ i = ... = u mb<i+ s = e F*. 

With the minimal encoder E mb which has v max memory 
elements, these 8+1 consecutive zeros of u mb (z) would result 
in the state vector <J mb ,t becoming zero for all time instants 
from i + u max to i + 6 + 1, i.e., 

Tmb,i+i/ max = &mbd+v max + l = ■■■ = CmM+(5+l =06 fl^n-"" 

With v max < 5, the path traced by Vu^g] traces at least one 
zero cycle on the trellis corresponding to the minimal encoder. 
This concludes the proof. ■ 
We shall now prove the bound on a. 
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Theorem 1: The slope a of a rate b/c convolutional code 
C with degree 5 is lower bounded as 



a > 



1 



1 



Proof: First we note the fact that every path in the state 
transition diagram is either a cycle or a part of a cycle. By 
Lemma [2] the path traced by any partial codeword sequence 
with 5 + 1 consecutive zero components in the state transition 
diagram of the minimal encoder would have a cycled around in 
the zero state at least once. The definition of a excludes a cycle 
around the zero state, and therefore paths (partial codeword 
sequences) which have 6 + 1 consecutive zero components 
cannot be considered to measure a since they would ultimately 
result in a zero cycle. However, Lemma|2]also implies that any 
path in the state transition diagram which does not include 
the zero cycle must therefore accumulate at least 1 Hamming 
weight in every 6 + 1 transitions. Thus we have proved that 
the slope is lower bounded as: 

1 

a > -. 

-6 + 1 



VII. Discussion 

The performance of CNECCs under the BSC edge error 
model has been analyzed using theoretical bounds and sim- 
ulations. A sufficient upper bound on the edge cross-over 
probability p e has been obtained, so that if p e is below this 
bound, the complexity of analysis can be reduced greatly 
by considering only single edge network-errors. Codes with 
better distance spectra and those with good slopes are seen 
to perform well under different conditions on the cross-over 
probability. A lower bound on the slope of any convolutional 
code is also obtained. Several interesting problems remain in 
this context including the following. 

• Studying the soft-decision decoding performance of 
CNECCs. 

• Constructions of convolutional codes with good slopes. 

• In large networks, error probabilities at the sinks could be 
large even for negligible p e values. It would be interesting 
to look at the existing network error correction schemes 
for such networks, and compare them with schemes 
which involve coding over smaller subnetworks. 
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Appendix A 
Convolutional codes-Basic Results 

We review the basic concepts related to convolutional codes, 
used extensively throughout the rest of the paper. For q, power 
of a prime, let F g denote the finite field with q elements, ¥ q [z] 
denote the ring of univariate polynomials in z with coefficients 
from Fg, ¥ q {z) denote the field of rational functions with 
variable z and coefficients from ¥ q and F g [[z]] denote the 
ring of formal power series with coefficients from F g . Every 
element of F g [[z]] of the form x(z) = X)i=o x i z% ^ x i ^q- 
Thus, ¥ q [z] C F 9 [[z]]. We denote the set of n-tuples over 
F 9 [[z]] as F"[[z]]. Also, a rational function x(z) = |||y with 
b(0) ^ is said to be realizable. A matrix populated entirely 
with realizable functions is called a realizable matrix. 

For a convolutional code, the information sequence u = 
[uo,Ui, ...,Ut] (ui £ ¥ q ) and the codeword sequence (output 
sequence) v = [vq, V\, v t ] (Vi £ F£) can be represented in 
terms of the delay parameter z as 



"0) = '■ 



and 



v(z) =^Vi 



i=0 



i=0 



Definition 4 ( M4V ): A convolutional code, C of rate 
b/c (b < c) is defined as 

C = {v(z)E¥< q [[z]] | v(z) = u(z)G(z)} 

where G(z) is a bx c generator matrix with entries from ¥ q (z) 
and rank b over ¥ q (z), and v(z) being the codeword sequence 
arising from the information sequence, u(z) € F^[[z]]. 

Two generator matrices are said to be equivalent if they 
encode the same convolutional code. A polynomial generator 
matrix l[l"4l for a convolutional code C is a generator matrix for 
C with all its entries from F g [z] . It is known that every convo- 
lutional code has a polynomial generator matrix [14|. Also, a 
generator matrix for a convolutional code is catastrophic lfl"4l 
if there exists an information sequence with infinitely many 
non-zero components, that results in a codeword with only 
finitely many non-zero components. 

For a polynomial generator matrix G(z), let gij(z) be the 
element of G(z) in the i th row and the j th column, and 

Vi := max deg(gtj(z)) 
j 

be the i th row degree of G(z). Let 

b 

8=1 

be the degree of G(z). 

Definition 5 ( H14V ): A polynomial generator matrix is 
called basic if it has a polynomial right inverse. It is called 
minimal if its degree S is minimum among all generator 
matrices of C. 

Forney in |16| showed that the ordered set {i^i, ^2, ^fe} 
of row degrees (indices) is the same for all minimal basic 
generator matrices of C (which are all equivalent to one 
another). Therefore the ordered row degrees and the degree 
S can be defined for a convolutional code C. Also, any 
minimal basic generator matrix for a convolutional code is 
non-catastrophic. 

Definition 6 ( H14\l ): A convolutional encoder is a physical 
realization of a generator matrix by a linear sequential circuit. 
Two encoders are said to be equivalent encoders if they encode 
the same code. A minimal encoder is an encoder with the 
minimal number of memory elements among all equivalent 
encoders. 

Given an encoder with 6' memory elements for the code 
C, we can associate a vector cr t G whose components 
indicate the states of the 6' memory elements at time instant 
t. 

The weight of a vector v(z) € F^[[z]] is the sum of the 
Hamming weights (over ¥ q ) of all its F^ -coefficients. Then 
we have the following definitions. 

Definition 7 ( H14V ): The free distance of a convolutional 
code C is given as 

dfree(C) = vain {wt(v(z))\v(z) G C, v(z) ^ 0} 



